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Abstract 

The quantum mechanical transition amplitudes are calculated perturbatively on the basis 
of the stochastic quantization method of Parisi and Wu. It is shown that the stochastic scheme 
reproduces the ordinary result for the amplitude and systematically incorporates higher-order 
effects, even at the lowest order. 
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1 Introduction 



Since its proposal, the stochastic quantization method of Parisi and Wu0 has been 
widely applied to various fields in physics 0, quite straightforwardly at relatively early 
stages of its development and with some conceptual and/or technical innovation at 
later times. A scheme for obtaining the quantum mechanical transition amplitude on 
the basis of the stochastic quantization (SQ) was proposed by H. Hiiffel and one of 
the present authors a few years ago|| and is considered to belong to the latter type 
of applications: One may easily realize the difficulty in obtaining an amplitude from a 
probability, on which the stochastic process, even if fictitious in this scheme, is essentially 
based. In this paper, this scheme|3|] is applied to a nonlinear system, for which no exact 
expression of the transition amplitude is known and we have to rely on a perturbative 
treatment. The purpose of the paper is to show how a perturbative treatment is possible 
within this scheme, to confirm its consistency with the ordinary one and then to find 
possible advantages over the latter. 

After brief reviews of SQ in Sec. |2] and of the stochastic scheme for the transition 
amplitude in || in Sec. we describe how to treat nonlinear systems within this scheme 
and propose a perturbative treatment in Sec. [|. Section |5| demonstrates some details 
of the calculation of the transition amplitudes for a quantum mechanical anharmonic 
oscillator. It will be clear that our treatment is equivalent to a perturbative expansion 
of the logarithm of the amplitude and therefore the amplitude thus obtained includes 
automatically a part of the higher-order contributions as an exponential form. The last 
section, Sec. ||, is devoted to a summary. 

2 Review of the Stochastic Quantization Method 

Let us briefly outline the ordinary prescription of the stochastic quantization method 
of Parisi and Wu[l|, ||]. This quantization method is so designed as to give quantum 
mechanics as the thermal-equilibrium limit of a hypothetical stochastic process. For 
this purpose, the dynamical variable x(t)f\ is assumed to be a stochastic variable x(t, s) 
with respect to a newly-introduced time, called the "fictitious time" s. Its dynamics is 
given by the Langevin equationP] 



where S is the classical action of the system and r] is the Gaussian white noise, whose 
statistical properties are characterized by 



tional of the noise rj, calculate the equal-time correlation function {x(t\, s)x(t2, s) ■ ■ ■) 

■"■For simplicity, we exclusively consider one-dimensional quantum mechanical systems here. 
2 We work in the Minkowski stochastic quantization schemeB throughout the paper. 




(2.1) 
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by means of ( |2.2| ), and take the equilibrium limit s — > oo. What we thereby obtain is 
shown to correspond to the ordinary vacuum expectation value (0\Tx(ti)x(t 2 ) • • • |0). 

This correspondence is most clearly seen in the Fokker-Planck picture. In this picture 
the stochastic average (• • ■) is expressed as a functional integral 

(x(t u s)x(t 2 , s) ■ ■ •) = J Vx xit^xitz) ■ ■ ■ P\x; s] (2.3) 

with a (quasi-)probability distribution functional P which obeys the Fokker-Planck 
equation 

d „ r r°° , 8 f 5 , 8S \ „ r 

—P[x; s] = / dt—— — — - i— - P[x; s] 2.4 
ds 1 1 J-oo 5x(t) \5x(t) 5x(t) J 

and is normalized as (1) = J VxP[x; s] = 1. Clearly, the stationary solution of this 
equation is given by e lS , which also serves as the equilibrium distribution if we adopt the 
so-called ie-prescription^J. Under this prescription, the equal-time correlation function 
(|2.3j ) approaches, as s — > oo, the Feynman path integral, which is nothing but the 
vacuum expectation value in the canonical (operator) formalism: 

Vx x(ti)x(i 2 ) • • • e iS 

lim (x(t 1: s)x(t 2 , s) ■ ■ •) = J - j = (0|Tx(*i)z(t 2 ) ■ ■ ■ |0). (2.5) 

Vxe lS 



3 The Stochastic Formula for Transition Amplitude 

The transition amplitude is one of the fundamental quantities in quantum mechanics, 
for the amplitude plays a central role there. On the other hand, as one readily notices, 
the above stochastic scheme apparently provides us with those expectation values that 
are normalized in the sense of (|2.5|) . This can be thought of as rooted in the essential 
property of the stochastic process, where the probability (density), not the amplitude, 
governs the dynamics. Therefore it does not seem trivial to derive amplitudes within 
the framework of the stochastic formalism of Parisi and Wu. 

Let us first try to calculate the correlation function, under the boundary conditions 

x(tp, s) — x-p, x(ti, s) = x\. (3.1) 

These conditions are considered to be one of the simplest extensions of those for the 
quantum mechanical transition amplitude (xp, ^fI^Ij ^i)- It would be natural to expect 
that the correlation function under the above boundary conditions (x(ti, s)x(t 2 , s) ■ ■ -)b c 
is expressed in the Fokker-Planck picture as 



rx(t F )=x F 

(x(t u s)x(t 2j s)---) hc = Vxx(t 1 )x(t 2 )---P[x;s] (3.2) 

Jx(ti)=xi 



with the normalization condition (l)bc = 1- The quantity with the subscript be is to be 
evaluated within SQ under the above boundary conditions ( |3.1| ). This normalization 
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condition implies that in the s 
following normalized quantity 



oo limit, the correlation function approaches the 



lim (x(ti, s)x(t 2 , s) ■ ■ -)bc 



x(t F )=x F 
x(ti)=xi 



T>X x{t\)x(t2) 



x(tp)=x F 
x(ti)=xi 



Vx e 



iS 



{x F , t F \Tx(ti)x(t 2 ) ■ ■ ■ \x h ti) 
(x F ,t F \x h ti) 



(3.3) 



It is such normalized expectation values ( |3.3| ) that we can directly calculate within 
this scheme. Remember that what we are interested in here is not the vacuum expec- 
tation value, but the transition amplitude (x F , t F \x\, t{) itself or, in other words, the 
normalization factor (i.e. the denominator in ( |3.3|) ) itself. 

In 0], a formula for the transition amplitude within the framework of SQ is proposed 
and applied to several solvable cases. The formula is based on the following relations 
in quantum mechanics 



and 



d 

-i-^-{x F ,t F \x h ti) = (x F ,t F \H(ti)\x h ti) 



d 

i-^-{x F ,t F \xi,ti) = (x F ,t F \p(ti)\xi,ti), 



(3.4a) 



(3.4b) 



where H(t\) and p{tj) are the Hamiltonian and the momentum operators, respectively. 
From flOaD we get 

— ln(x F , t F \x h ti) = i ; — : -r = i(H(ti)) hc , (3.5) 

Oti \^F; t F \Xi, ti) 



that has the solution 



(x F , t F \xi, ti) = cexp 



i p dti(H(ti)) 



be 



(3.6a) 



with a factor c dependent on x F , xi and t F . Along the same line of thought, we obtain 
from (|3~4bD 



(x F , t F \xi, ti) = cexp 



J,'! 



i / dxi(p(ti)) 



be 



(3.6b) 



with a factor c dependent on x F , t F and tj. The transition amplitude (x F ,t F \xi,ti) has 
thus been related to the (normalized) expectation value of the Hamiltonian (H(ti))h c 
and to that of the momentum (p(ti))bc, which are obtainable as the equilibrium limits 
of (H(t h s)) hc and (p(t h s)) hc in SQ.0 

3 Equation ( |3.6a ) is the formula presented in j3|. On the other hand, though it is not explicitly 
presented there, ( |3.6b| ) is surely used in the practical calculations. 
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Notice that these relations are nothing but two of the expressions of the ordinary 
variational principle and there are two more formulae, i.e., the formulae relating the 
transition amplitude to (-ff(tp))bc an d {p(ip))bc- For most systems, however, the above 
two formulae (|3.6| ) are sufficient to determine c and c except for a constant factor inde- 
pendent of tp, ti, xf and x\. In fact, we can determine the xi-dependence of c from ( |3.6b| ) 
and the tp-dependence of c from ( |3.6a| ) . Furthermore, their xp- and tp-dependences are 
fixed from the fact that (xp, £f|xt, tj) is a function of T = ip — ti (translational in- 
variance) and is a symmetric function of xp and x\ (time- reversal invariance). The 
constants c and c are thus completely determined if we fix the remaining factor by 
requiring (xp, tp\xi, h) to approach a Dirac (^-function <5(xf — x\) as T — > 0. 



4 A Perturbative Treatment 

The above formulae (|3.6|) have been applied in [[| to linear systems to derive transition 
amplitudes within the framework of SQ. It is demonstrated that the stochastic scheme 
can reproduce the correct results for the amplitudes. In the following, we shall develop 
a perturbative method in order to treat nonlinear systems within this scheme. 

Since all that is needed in deriving the amplitude in SQ are the expectation values 
of H and p, according to ( |3.6a|) and (|3.6b|) , a natural strategy would be to evaluate 



these quantities as a power series in the coupling constant. For systems in which the 
Hamiltonian H is composed of a kinetic part p 2 /2M and a potential part V(x), what 
essentially remains to be evaluated is the expectation value of the former, owing to the 
boundary conditions (|3.1|) 

(i/(tl))bc = ^:(p 2 (tl))bc + V(x I ). (4.1) 

We have to perturbatively evaluate the expectation values of momentum and squared 
momentum. Though this can be done within the framework of the phase space formal- 
ism of SQ[||, |3|, we prefer to work in ordinary configuration space for simplicity. It is 
not difficult to derive the following relations 

(p(*i))bc = M lim dt lim (x(t, s)) hc (4.2a) 

t — >t\ s — >oo 



and 



(p 2 (ti))bc = M 2 lim d tl d t2 lim (x{h, s)x(t 2 , s)) bc , (4.2b) 

tl,r2 — >ti s — *oc 



on the basis of the phase-space Langevin equations.^ Thus the problem is reduced to 
the evaluation of the correlation functions, (x(t,s)) bc and (x(ti, s)x(t 2 , s))bc : We solve 
the Langevin equation ( |2.1| ) perturbatively and get its solution x(t, s) as a power series 
in the coupling constant. As in the solvable cases treated in ||, the boundary conditions 

4 The second relation just corresponds to Feynman's time-splitting procedure in evaluating the 
average square velocity ||. 
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( |3.1| ) shall be taken into account by the classical solution and the remaining part (to 
be quantized in terms of the Langevin equation) is subject to the "zero boundary 
conditions" and is expressed as a Fourier series. 

It is worth mentioning here that the perturbative treatment on the basis of the 
formulae p.6a| ) and ( |3.6b| ) differs from the ordinary one: What is really evaluated 
perturbatively is not the amplitude itself, as in the ordinary case, but its logarithm. 
Therefore we can expect that even the lowest-order value of (H(ti)) hc or (p(ti)) hc will 
systematically yield higher-order contributions to the amplitude. 

5 An Example: Anharmonic Oscillator 

Let us consider an anharmonic oscillator described by the action 



S = J dt [t,Mx 2 - -Mu) 2 x 2 - |a; 4 J (5.1) 

and calculate the transition amplitude (x F , t F \xi, ti) according to the above stochastic 
scheme. The Langevin equation ( |2.1| ) reads 

d s x{t, s) = -iM(d 2 + u 2 )x(t, s) - igx 3 (t, s) + rj(t, s) (5.2) 

and we write its solution as the sum of the solution for a harmonic oscillator (g = 0) 
x^ and a remaining part 

x(t,s) = x {0 \t,s) + x'(t,s). (5.3) 
The free solution x^ is composed of two terms 

x^(t,s) = x^\t)+x^(t,s), (5.4) 

where x^ is the solution of the classical equation of motion 

(d 2 + oo 2 )xf(t) = with x^\t F ) = x F , x^(h) = xi (5.5) 
and x'q that of a free Langevin equation 

d s x { $ (t, s) = -%M{d 2 t + uj 2 )x { q ] (t, s) + 7](t, s) with x^ (t F ) = x§ ] (h) = 0, (5.6) 

where the last conditions are called "zero boundary conditions." They are readily cal- 
culated to be | 

2vi (t) = = \x F sin u(t — ti) + X\ sin u(t F — t)\ (5.7) 

smuT 

and 

4° } (M) = f tF dt' r ds'G { V(t,t f ;s-s'W,s'), (5.8) 

Jt\ J —oo 
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Fig. 1. Graphical representation of the perturbation series of x. 



where 



Gg (t, t'; s -s') = e(s-s')^± sin ^ (t - h) sin ^ (f - h)e 

1 71=1 



(W) 



(5.9) 



is the free Green function which respects the "zero boundary conditions" The 
retarded Green function G b ^ satisfies 

d s + iM(d* + u 2 )] G^{t,t'- s-s') = 6(t - t')8(s - s') (5.10a) 
and the boundary conditions 

G$(t, t'- 8-s') = for t (or t') = t F , tj or s < s'. (5.10b) 



Here we have set the initial fictitious time at so = — oo, in order to achieve the long-time 
limit automatically. This choice makes any initial- value dependence of the solution Xq* 
irrelevant, for the system is considered to be in equilibrium already at finite s. Now we 
solve the Langevin equation Q5.2|) to get a recursive relation for x(t, s) 



x(t, s) = x§\t) + x§\t, s)-ig dt' / ds'G®(t,t'; s - s')x 3 {t', s 1 ), (5.11) 

JU J —oo 



from which we obtain the perturbation series 

r _ T (0) , _(0) _ f tF df > [°° H j a (0) (AO) 3 , o T (0) 2 (0) (0) (O) 2 (0) ; 

x — x cl X Q l y / ut / uil ^bc I x cl ox cl X Q ox cl X Q < X Q 

Jtj J —oo \ 



+ 0(g 2 ). 
(5-12) 

We represent the series ( |5.12| ) diagrammatically as in Fig. [TJ, where x^\ Gjjj, an< ^ ^ 
are denoted by a line with a dot, a line with an arrow, and a cross, respectively, and at 
each vertex, a factor —ig and integration variables t' and s' are attached. The stochastic 
average over rj (|2.2|) combines all the crosses in pairs in all possible ways. We also have 
to integrate over the vertex times t' and s' (see (|5.12 )). 

According to the formula (|4.2a|) , we need to evaluate (x)b c in order to get (p(ti))b c > 
to which the three diagrams in Fig. ^| contribute up to O(g), where a simple line without 
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• 



Fig. 2. Stochastic diagrams contributing to (x)b c up to 0(g). 

any arrow or dot denotes the free two-point correlation function 

D^(tAs-s>) 
= (x^(t,s)x^(t',s')) hc 

= 2 dt x r ds 1 G^(t,t 1 ;s-s 1 )G^(t',t 1 ;s , -s 1 ) 



oo 

2 



2i 1 riir nn -iMu) 2 

- - nrTn > sin — r(t — h) sin — -It — ti)e 

Mu 2 T ^ _ x T ' T ' 



l*-s'| 



(5.13) 

These three diagrams correspond to the following terms in {x)^ c 
(x(t,s)) hc 

= 4\t) ~ W f F dt' r ds'G^(t, t>- s - s') (xf\t>) + 3xf(t')D^(t', t'- 0)) 

Jt\ J —oo \ / 

+ 0(g 2 ). (5.14) 

A relation between G^} and (see the Appendix) enables us to perform the integra- 
tion over s' and yields 

(x(t)) bc = x^(t) - igfj dt'A hc (t,t') (x { f(t>) + 3x < S\t')A bc (t', f )) + 0(g 2 ), (5.15) 

where A bc (t,t') = D^(t,t'; 0) and its explicit form is|| 

Ahc(t,t') = — —smu(tF — max(t,t'))smu(mm(t,t') —ti). (5.16) 

Muj sm ujT 

On the other hand, the expectation value (p 2 (ti))b c is calculated from the correlation 
function (xx)b c through the relation ( |4.2b|) . At the lowest order, there are three types 



of connected diagrams, shown in Fig. §, contributing to (xx) bc and we have, after the 
integration over s', 

(x(t)x(t')) hc = (x(t)) hc (x(t')) hc + A hc (t,t') 
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Fig. 3. Connected stochastic diagrams contributing to (xx)bc up to 0(g). 



- 3ig dhA hc (t, h) ^ )2 (tx) + A bc (t!, tjj) A hc (t u t') 

+ 0(g 2 ). (5.17) 

We are now in a position to write the explicit form of the transition amplitude 
(xp, tp\xi, ti) for the present system. Since the Hamiltonian of this system is given by 



H 



1 



we have 



(H(ti)) 



2M 
1 



p z + -Mu 2 x 2 + -x 
2 4 



be 



(p 2 (ti))bc + -Mu?x\ + g ~x\, 



(5.1* 



(5.19) 



2M U " v " i//u " ' 2 1 ' 4" 
where the last two terms are a consequence of the boundary conditions ( |3.1| ) for x. 
From (|57T5|) and fl5T/D we calculate (p(ti)) bc and (H(ti)) hc through (]57T9|) and (g3|), 
and insert them into the formulae (^.6|). In this way, we obtain (xp, ^f|^i, *i) m the 
exponential form 



where 

and 
Si = 



x F ,t F \x h ti) 
Muj 



Muj 



2 sin tuT 



2m sin ujT 

„2 , 2 



3 iS +i5i+O(g 2 ) 



x ¥ + x Y ) cos ujT — 2xpX\ 



(5.20a) 
(5.20b) 



9 



32M 2 uj 3 sin 2 cjT 
3 

M ^16Mw 2 sin 3 uoT 



1 



5 



32co>sin ujT 



{?>ujT — 3 sin ujT cos cuT — 2c<jT sin 2 ujT 

[/ Xp + x 2 ) \3ujT cos luT — 3 sin co>T + sin 3 tuT 

— 2xpX\ [?>ujT — 3 sin uT cos uT — 2u>T sin 2 ujT 
x F + xf) (SluT — 3 sin uT cos ujT — 2 sin 3 cuT cos uoT 



— A(x ¥ xi + x F x\) (3ujT cos tuT — 3 sin ujT + sin 3 ujT 



+ 6xpa; 2 ( 3a; T — 3 sin ujT cos ujT — 2cuT sin 2 u;T 



(5.20c) 
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Notice again that the constant c in ( |3.6a|) contains enough information to fix the 
T — tp — t\ dependence of c in ( |3.6b ). On the other hand, the latter can fix the 
xp(x\) dependence of the former. The remaining factor has been fixed from the normal- 
ization condition at tp — ti- The result ( |5.20| ) coincides with the ordinary perturbative 
expression for the transition amplitude if it is expanded as a power series in the coupling 
constant g, as it should be. 



6 Summary 

In this paper, we have proposed a perturbative treatment of transition amplitudes within 
the stochastic quantization scheme and have applied it to an anharmonic oscillator. Our 
treatment is based on the Langevin equation: We solve it iteratively to get its solution x 
perturbatively under given boundary conditions. This enables us to calculate (H(ti))b c 
and (p(ti))bc as power series in g, from which the transition amplitude is obtained by 
making use of the formulae (|3.6| ). Even though up to the order of interest (0(g), in 
this case), they coincides with each other, the expression (|5.20| ) is different from the 



one usually obtained, for example, in the path-integral formulation: In this scheme, it 
is the argument of an exponential function that is expanded as a power series in g. In 
this sense, we can say that the present stochastic scheme systematically incorporates 
higher-order effects. 

It is worth stressing again that the quantity to be calculated perturbatively in this 
scheme is not the amplitude itself, but its logarithm. This is rooted in the exponential 
form of the formulae ( |3.6| ). These formulae are quite remarkable, in the sense that all 
quantum fluctuations are put together in one exponent (/ (H(t\))^ c dti or / (p(t\)\ c dxi) , 
which is contrasted to the usual situation in the path-integral quantization where fluc- 
tuations manifest themselves as different paths. Though it is not easy to get an intuitive 
physical image (if any) of these exponents, it might be helpful to remind that a similar 
relation exists between the ordinary generating functional and its counterpart for con- 
nected diagrams. The study, along such a line of thought, might lead us to a deeper 
understanding of the stochastic scheme itself and shed new light on the meaning of 
quantization. 
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Appendix 



We shall prove the relation 

ds'G i V(t,t';s-s') = A hc (t,t'), 

J — oo 

by means of which we get (|5.15|) . Here A bc is defined as 

A hc (t,t')= D^(t,t';s-s') 

with 

D^(t, t>- s - s') = (x ( °\t, s)x$\t', S '))bc 
This function has the following explicit form (see (|2~2"|), (|5l8| ) and (|5T9| 



(A.l) 

(A.2) 
(A.3) 



2i 



Mu 2 T 



1 717T , WK . -iMui 2 

— sin — (t — ti) sin — (£ — ti)e 



It is easy to recognize the existence of a relation between in flOl) and D^J in (|A.4| ) 



)(0) 



.(A.4) 



G^(t,t'; s-s') = 9{s - s')^D^{tA s - s') 



<9 



)(0), 



(A.5) 



Hence the integration in ( |A.1| ) becomes trivial. The above relation is nothing but 
the "golden rule," as first observed in and later interpreted as a realization of the 
fluctuation dissipation theorem[§]. 
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